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ABSTRACT 

Multiple image gravitational lens systems, and especially quads are invaluable in deter¬ 
mining the amount and distribution of mass in galaxies. This is usually done by mass 
modeling using parametric or free-form methods. An alternative way of extracting 
information about lens mass distribution is to use lensing degeneracies and invariants. 
Where applicable, they allow one to make conclusions about whole classes of lenses 
without model fitting. Here, we use approximate, but observationally useful invariants 
formed by the three relative polar angles of quad images around the lens center to show 
that many smooth elliptical-l-shear lenses can reproduce the same set of quad image 
angles within observational error. This result allows us to show in a model-free way 
what the general class of smooth elliptical-|-shear lenses looks like in the three dimen¬ 
sional (3D) space of image relative angles, and that this distribution does not match 
that of the observed quads. We conclude that, even though smooth elliptical-l-shear 
lenses can reproduce individual quads, they cannot reproduce the quad population. 
What is likely needed is substructure, with clump masses larger than those responsible 
for flux ratio anomalies in quads, or luminous or dark nearby perturber galaxies. 

Key words: gravitational lensing: strong - galaxies: fundamental parameters - dark 
matter 


1 INTRODUCTION 

The field of multiple image gravitational lensing was born 
about 35 years ago with the discovery of Q0957-I-561 ( |Walsh| 
|et al.|1979 1 and PG1115-I-08 ( Walsh et al. [198^ 1, which were 
followed by discoveries of many more doubly and quadruply 
imaged quasars. The theoretical understanding of these sys¬ 
tems grew up alongside the growing body of observations. 
The earliest works found that simple parametric forms for 
the mass distribution in the lens can account for the systems’ 
observables ( Young et al.|1980 l, and that image properties, 
namely positions, shapes, time delays and flux ratios can be 


derived from the suitably defined lensing potential (Bland- 
ford fc Narayan||1986 Schneider et al. 19921. 


It was also soon realized that strong lensing has degen¬ 
eracies and invariants, and that these theoretical insights 
prove useful when lensing is used as a tool. For example, 
mass sheet degeneracy, and some other degeneracies where 
different lenses reproduce exactly the same image observ¬ 
ables have been well studied (Falco et al.||1985 Gorenstein 


et al.|1988 Saha|2000 Liesenborgs & De Rijcke 20121, espe¬ 
cially in relation to mass modeling, where degeneracies affect 


the conclusions about the mass profiles of lenses. The first 
lensing invariant discovered was the magnification invariant. 
It states that the sum of signed magnifications for all lensed 
images of any source within the caustic of a given lens is a 
constant. This constant can be the same for more than one 
lensing potential, and is often independent of the parameters 
characterizing the pote ntial (|Witt fc Mao' 1995 Dalal|2005 1. 
I Witt fc Mao| ( [2000| ) and |Hunter fc Evans||200T| showed that 
image positions can also be used in forming invariants. Aside 
from their intrinsic theoretical value, the invariants allow us 
to rule out some lens models without even fitting the lens. 
To quote Hunter & Evans (20011, “The major application 


of lensing invariants is to shortcut the modeling process.” 

Because observational uncertainties on the image ob¬ 
servables are never zero, not only exact but also approximate 
degeneracies and invariants are important for the practical 
work. For example, [Schneider fc Slu'se] ( |2014| start with an 
exact invariance transformation that applies to axisymmet- 
ric lenses, and show that it survives in an approximate, but 
still useful form when applied to a wider range of lens mod¬ 
els. 


* E-mail: woldesenbet@physics.umn.edu (AGW); 

llrw@astro.umn.edu (LLRW) 


Our work in the present paper is in the same vein: we 
study near invariants that provide useful insights without 
mass modeling, and have important consequences for mass 
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substructure in lenses. Our analysis is of quads; it involves 
image positions only, and does not rely on magnifications. 

The relative image locations of a quad are specified 
by six numbers, say, (x, y) coordinates of any of the three 
images with respect to the 4th. Lens mass modeling— 
parametric or free-form—attempts to reproduce all six num¬ 
bers at once. In Williams et al. (|2008 1 and Woldesenbet 


& Williams (2012[) (hereafter WW12) we started develop¬ 


ing a new way of looking at quads, that considers only the 
three polar angular coordinates of quads, and disregards the 
three r coordinates. We showed that information about the 
mass distribution of certain lens models contained in the az¬ 
imuthal, or angular coordinates of images around the center 
of the lens is almost independent of that contained in the ra¬ 
dial coordinates, and that the angular coordinates of images 
show approximate degeneracies. Specifically, if a lens with 
double mirror symmetry can reproduce the relative images 
angles of a quad, then all lenses with double mirror sym¬ 
metry can reproduce the relative image angles of that quad, 
though for different locations of the (unobservable) source. 
Here we continue developing the theory of relative image an¬ 
gles of quads by extending it to lenses represented by purely 
or approximately elliptical mass profiles with external shear. 
The lenses in this general class of mass models are impor¬ 
tant to study because they fit observed quads well, and are 
common in parametric lens modeling. 

Since we study a wide class of such lenses, and identify 
similarities, we are able to make conclusions that are inde¬ 
pendent of a specific lens model. Our analysis suggests that 
even though these common models are able to reproduce 
observed quads one at a time, they are unable to reproduce 
them as a population. Instead, substructured lenses or lumi¬ 
nous or dark nearby perturbers are likely needed. 

Because of their importance in the ACDM cosmolog¬ 


ical model (Moore et al. 1999 Klypin et al. 19991, sub¬ 


structure has been the focus of many recent papers. The 
main method employed for finding substructure uses flrrx 
ratios of close pairs or triples of images of quad lenses. In 
the absence of substructure clumps, i.e. when the potential 
is smooth, the magnifications of these images obey certain 


relations (Blandford & Narayan 1986 Schneider & Weiss 


1992 ). Substructure induces deviations from these relations, 
or anomalies ( Mao fc Schneider|1998 Metcalf fc Zhao|2002 1 . 
The substructure finding methods based on flux anomalies 
and image positions are complementary because the former 
are sensitive to small clumps located close to the images, 
while the latter are less easily perturbed and respond to 
larger clumps, > 10 ® — IO^^Mq, located anywhere around 
the Einstein radius. The two methods also differ in other re¬ 
spects. Image fluxes are subject to different effects, depend¬ 
ing on wavelength, like microlensing by stars in the optical 
and near infra-red, and propagation effects in the radio, and 
it is still unclear how major a role these factors play in the 


observed flux anomalies (Xu et al. 20151. Image positions 
are affected only by relatively massive clumps. 

This paper is organized as follows. In Section we de¬ 
fine our classification of lens types based on the caustics they 
produce. In Sections and we discuss the detailed prop¬ 
erties of lenses that belong to what we call Type II lenses, 
which are the focus of this paper. Readers interested primar¬ 
ily in substructure can skip directly to Section Section 
summarizes our findings. 


2 TYPE II LENS MODELS 

Parametric modeling, where parameters such as ellipticity, 
shear, position angle, etc. identify a lens model, is a com¬ 
monly used method to theoretically represent gravitational 
lenses. Simple models obeying certain symmetries are used 
to fit observed properties of lens systems, such as image 
positions, time delays, and magnifications. The symmetries 
could be either in the isopotential or isodensity contours of 
lenses with different radial density profiles. Examples of such 
axisymmetric models that are widely used are Elliptic Mass 


Distribution (EMD) and Elliptic Potential (EP)(Kassiola & 
Kovner 19931 with different radial profile such as Isother¬ 


mal and NEW. The amount of constraining information that 
can be obtained from a given set of observed images is very 
limited, and inadequate to fully describe the details of the 
lensing mass. This means the lens equation is severely un¬ 
der constrained allowing it to admit multiple solutions. But 
any symmetry in how images are distributed on the lens 
plane is a direct consequence of the underlying symmetry 
of the lens and the source position in the source plane. For 
example, any axially symmetric lens, independent of the ra¬ 
dial density profile, would give rise to an Einstein ring for 
a central source. Therefore categorizing and studying lenses 
that share sets of symmetries using parameters is a practical 
approach that does not contradict observational data. 

In previous work by the authors (WW12) three classes 
of lenses were introduced. The classification is based on az¬ 
imuthal symmetries of isopotential or isodensity contours of 
lenses, and is independent of the radial density profile. Lens 
models obeying twofold symmetry (symmetric about two or¬ 
thogonal axes) in the lens plane were termed Type I lenses, 
and breaking this symmetry once resulted in Type II lenses. 
Type III encompassed all other models. 

A more precise categorization, which we will adopt in 
this paper, can be achieved by using the symmetries of the 
diamond caustics. Type I lenses are those that give rise to 
caustics with twofold rotational symmetrj 0 and double mir¬ 
ror symmetries, while Type II lenses are the ones that give 
rise to caustics obeying only twofold rotational symmetry. 
In this new definition of Type I and II we also limit our¬ 
selves to lenses with diamond caustics and exclude those 
that produce caustics with higher order catastrophes. Such 
constraint is justified for our analysis which aims to under¬ 
stand the population of observed lens systems. Higher order 
catastrophes can produce lenses with more than five images, 
which are not observed in galaxy lenses. Type HI lenses in¬ 
clude everything except Type I and Type H, for example, 
substructured lenses. This paper examines Type H models 
in detail. 

Elliptical potentials or mass distributions of arbitrary 
ellipticity and density profile have twofold and double mir¬ 
ror symmetries. By introducing external shear at a non-zero 
angle with respect to the ellipticity axes, the double mirror 
symmetry is broken giving rise to Type H lenses. In this pa¬ 
per we look at Type Hs with power law and NEW density 
profiles. We identify a lens by three parameters; ellipticity, 
e, external shear, 7 , and the angle, /3, between shear and 


^ Objects with n-fold rotational symmetry look the same when 
rotated by =11^ degrees. Twofold symmetry is identical to inver¬ 
sion symmetry through the origin. 
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the ellipticity principal axis. By definition, /3 lies between 0° 
and 90°. 

Our basic analysis method uses the relative angular dis¬ 
tribution of the four images of quad^ about the center of 
lens. It is applicable to any lens system, such as quasar- 
galaxy or galaxy-galaxy, where the center of the lens is 
known and the four images are point like. This approach was 
introduced in Williams et al. (20081 and refined in WW12. 
The method is statistical in nature in the sense that it works 
with quads as a class to draw conclusions about the mass 
distribution in galaxy hosts of quad lenses. We will show 
that properties related to quad image angles are approxi¬ 
mately independent of the radial mass density profile of the 
lens, and that Type I and II lenses cannot account for the 
observed population of quads. 


3 POWER LAW POTENTIALS AS EXAMPLES 
OF TYPE II 


Type II lenses are a wide class. Our strategy is to study a 
few representative examples of Type IIs in detail, note their 
similarities and differences, and draw conclusions based on 
these. Our examples are chosen to be representative of the 
models used in the lensing literature and to resemble the 
real galaxy lenses. Furthermore, for computational ease, we 
use models whose lensing potentials are expressed by simple 
analytical functions. 

Power law lensing potentials are one example that sat¬ 
isfy all the above criteria. They have a general form, 

(l>{r, 9) = br°‘^/T^^ems{26) +cos(2[9 — P]). (1) 


In the above expression b is the normalization factor, e is 
an ellipticity parameter (henceforth ellipticity|^ and r is 
the sky-projected distance from the lens center. The second 
term is external shear of strength 7 oriented at an angle of 
P relative to the ellipticity’s principal axis. The exponent 
a in the first term is physically constrained to be between 
zero (q < 0 implies decreasing total enclosed mass with 
increasing r), and two (a > 2 implies that the density is 
increasing with r), but observations constrain it even further 
to be around 1 (Koopmans et al. 20061. Theoretically, we can 
use any values of 7 as long as the lens gives rise to quads, 
whereas e is limited to be between 0 and 1 by definition. 
In addition, e and 7 are also restricted from above by the 
requirement that the corresponding isodensity contours are 
not peanut shaped, which happens, for example, for a « 1.2 
and e « 0.25, or a « 1 and e « 0.15. 


3.1 Singular Isothermal Elliptical Potential with 
External Shear (SIEP+shear) 

In this section we use singular isothermal elliptical potential 
with external shear to demonstrate detailed properties of 


^ Technically, quads are five image configuration with one of the 
images invisibly sitting near the very bright center of lenses. For 
all purposes of this paper, the fifth central image is ignored. 

® Note that e is related to the standard definition of ellipticity 

(1 - axes ratio) as 1—^ = 1 — 



Type II lenses. SIEP-(-shear, given by eq. [^with a = 1, is 
chosen because it is the simplest form of power law potential 
that is semi-analytically treatable. Different combinations of 
e, 7 and /3 give rise to different diamond caustics (based on 
orientation, elongation of cusps and curving of folds) but all 
obey the twofold rotational symmetry that identifies Type 
II lenses. 

One interesting property of Type IIs concerns the isopo¬ 
tential contours. The main galaxy, i.e. first term in equation 
[^with a = 1 gives rise to contours all of which have the same 
ellipticity axis, but the addition of external shear results in 
slightly twisting isopotential contours, i.e. angular orienta¬ 
tion of the principal axis of each contour of the total lens 
potential changes with radius. This is true of the contours 
in the first panel of Figure]^ but is hard to see because the 
degree of twisting is small. 

We parametrize a given diamond caustic by angles 
formed by its diagonals and the ratio of the length of the di¬ 
agonals. We do not consider the curvatures of the folds, i.e. 
there could be two diamond caustics of the same diagonal 
ratios and angles but with different fold curvatures. In the 
case of Type I lenses, the caustic diagonals are perpendic¬ 
ular to each other while for Type II caustics, the diagonals 
form angles other than 90°; see the middle panel of Figurej^ 

Another interesting property of SIEP-(-shear is that the 
singularity at the center is not a true critical point; mapping 
it to the source position using the lens equation gives rise 
to a set of points forming an oval which acts like an oval 
caustic, and is called a ‘cut’ (Kovner 1987(. Even though 
this set of points does not satisfy the common definition of 
a caustic, where the determinant of the Jacobian of the lens 
mapping is zero, it still acts like one; crossing the cut results 
in a change of image multiplicity. 

Quads are formed when a source is within the diamond 
caustic and the cut. The angular distribution of the four im¬ 
ages about the center of the lens can be uniquely represented 
by three relative angles 9ij (the acute angle between image 
i and image j), where image i is the arriving image. Im¬ 
age ordering for synthetic lenses is always known. Images in 
observed quads can be correctly time ordered in most cases 
using image morphology ( Saha fc Williams|2003 (. 

In this paper we use the same set of relative angles as 
in WW12: ^ 12 , ^ 34 , and 623 . A 3D space can now be formed 
using these angles, where a point in the space represents 
a single quad. Given a lens characterized by (e, 7 ,/ 3 ), one 
can generate a large number of quads arising from differ¬ 
ent source positions. Numerical calculation shows that the 
normalization factor, 6 , has no impact on the angular distri¬ 
bution of the images. This also holds true for the NFW lens 
discussed later in this paper. The distribution of the corre¬ 
sponding points in the 3D angle space can be used as an 
alternative way of characterizing the lens. Quads produced 
by Type I lenses lie on a slightly curved surface, which we 
called the Fundamental Surface of Quads (FSQ) with a peak 
at ( 012 , 034 , 023) = (180°, 180°, 90°). Quads from all Type 
I lenses lie very close to the FSQ, making it a nearly invari¬ 
ant surface, and a useful reference. 

On the other hand, quads arising from Type II lenses 
form two separate surfaces, each usually confined to two 
different portions in the 3D angle space separated by the 
FSQ. Based on the positions relative to the FSQ, the two 
surfaces are identified as upper or lower surface. Like the 
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FSQ, each of these surfaces has two well defined edges which 
corresponds to 612 = 180° and 634, = 180°. These two edges 
of each of the two surfaces meet to form two peaks, one above 
and one below the FSQ peak as shown in the left panel of 
Figure Each surface sits at the farthest point relative to 
the FSQ at its peak, i.e. near 6^23 ~ 90°. The edges of the 
surfaces come closer to each other at lower values of angle 
023 while their middle parts stay farther apart (right panel 
of Figure [^. The two 623 angles of the peaks of the two 
surfaces of an SIEP+shear lens (or any Type II lens) are 
supplementary angles (i.e. they add to 180°). We use the 
023 value of the peak quad of the upper of the two surfaces, 
023 ,p, as an index to parametrize the surfaces of a given lens 
in the 3D space. 

3.1.1 Potential, caustic and 3D angles space 

As described above. Type II lenses can be characterized in 
three different spaces (Figure[^. The potential space, where 
e, 7 and /3 parametrize the lens, the caustic space where 
the angles between caustic diagonals and the ratio between 
their lengths characterize the caustic, and the 3D angles 
space where 023 ,p value of the peak of the top surface is the 
characteristic parameter. Now, we describe how the three 
spaces and their corresponding parameters are related. 

The central source. The source at the center of the 
lens corresponds to the center of the diamond caustic and is 
mapped to two pairs of images with each image of a given 
pair having the same arrival time. We arbitrarily choose one 
of the first arrivals as image 1 and the other image 2. The 
same is done with images 3 and 4. The two images with the 
same arrival time sit at 180° of each other on the lens plane, 
therefore 0 i 2 = 034 = 180°. For Type I’s, the line connecting 
opposite images (those with the same arrival time) are or¬ 
thogonal, so 023 = 90°, but this is not the case for Type IIs, 
where 023 is either acute or obtuse, depending on the choice 
of ordering, but the two possibilities are always supplemen¬ 
tary angles. The maximum angular separation of images 2 
and 3 is attained when the source is at the center of the caus¬ 
tic, and this separation decreases as the source moves radi¬ 
ally out from the center. In the 3D angles space, the central 
source corresponds to two degenerate quads, each located at 
the peak of each of the two surfaces, at (180°, 180°, @ 23 ,p), 
and (180°, 180°, 180° - 023 ,p). 

SIEP-l-shear has an interesting property connecting the 
caustic and the images arising from the central source. The 
two supplementary angles formed by the diagonals of a dia¬ 
mond caustic are same as the 023,p angles of the 3D space. 
The symbolic expressions for the angles formed by the diago¬ 
nals of the diamond caustics and 023 ,p are rather large when 
expressed in terms of the parameters e, 7 , and /3. There¬ 
fore, we compared them by calculating the difference be¬ 
tween them using Wolfram Mathematica for 10,000 differ¬ 
ent combinations of (e, 7 , and /?) which resulted in exactly 
zero up to 13 decimal place precision. (In addition, the ratio 
of the caustic diagonal lengths is the same as the ratio of 
the radial positions of images 2 and 3 from the center of the 
lens, but we do not consider the image distance ratios in this 
paper, so we will not explore this relation any further.) 

Sources on the caustic diagonals. The images of 
sources on the diagonals of the caustic form the two outer 
edges of the sheets, (180°, ^34, ^23) and {0i2, 180°, @23). Just 


like the central source, sources on the diagonals results in 
degenerate images in terms of arrival time. It is important 
to note that this degeneracy is only for source positions 
along the diagonals. All other sources within the caustic are 
mapped to quads that form the body of the 3D surfaces, 
with each image arriving at distinctly different time from 
the rest. Therefore, the two sheets are completely different 
and not a result of ordering choice of the equally arriving 
images. 

The rest of the sources. The rest of the two sheets 
in the 3D angles space comprises quads arising from sources 
within the body of the diamond caustic. As shown in the top 
panels of Figure there is a bifurcation about the diagonal 
of the caustic. Continuously crossing a diagonal, i.e. moving 
from one quadrant of the diamond caustics to another (for 
example, from the black region to the yellow, lighter shade, 
one) results in a jump between the two sheets of the 3D 
angles space. Sources on straight lines of constant source 
position angle in the source plane (various colored lines in 
the caustic shown in the lower left panel of Figure]^ result in 
quads that form non-crossing monotonic, but not necessarily 
straight curves on the 3D surfaces (lower right panel). 

The inversion symmetry of Type II lenses. The 
twofold (inversion) symmetry of the Type II lens diamond 
caustic also applies to the potential, and is further reflected 
in the 3D angles space. Source positions from the two oppo¬ 
site quadrants of the caustic generate quads that form one 
of the sheets. The remaining two quadrants produce the sec¬ 
ond sheet. In Figurej^the two yellow, lighter shade, (black) 
quadrants of the caustic shown in the upper left panel gener¬ 
ate the yellow, lighter shade, (black) surface in the 3D angles 
space shown in the upper right panel. So the inversion sym¬ 
metry reduces the number of distinct caustic quadrants from 
four to two. The existence of the two surfaces in the 3D space 
of relative image angles is the most important characteristic 
of Type II lenses. 

Note that the properties described in this subsection 
are not unique to SIEP-|-shear lenses, but are common to 
all Type IIs. The inversion symmetry of Type II, which is 
their defining property, naturally predicts the two surfaces 
in the 3D angles space. We have tested the existence of the 
two sheets on many more lenses than are presented in this 
paper (for example, several Sersic type models), and have 
confirmed the qualitative behavior of the two sheets. At large 
023 the two surfaces are always found on the opposite sides 
of the FSQ in the 3D angles space (left panel of Figure]^, 
and the separation of the two peaks depends on how much 
the lens potential deviates from that of Type I. At small 
023 the edges of the two sheets approach each other and the 
FSQ (right panel of Figure [^. 

3.1.2 Near degeneracies of lens models with different sets 
of (e, 7, d) 

As mentioned earlier, we chose to parametrize the distri¬ 
bution of quads in the 3D angles space by 023,p, the angle 
between the 2nd and 3rd arriving images of the source lo¬ 
cated at the center of the lens. The full distribution of quads 
from a single lens are the two surfaces or sheets (right panel 
of Figure]^ or Figure]^. One may ask if it is adequate to 
represent the surfaces with just a single point, 023 ,pi 

Figure shows contour surfaces of constant 023,p in 7 
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vs. e vs. P space. Any two points on a given surface represent 
two different lenses in the potential space, each character¬ 
ized by a different set of (e, 7 , P), but sharing the same 
^ 23 ,p. The two contour surfaces shown have 023,p = 91° and 
023,p = 92°. In the caustic space these lenses have some 
commonalities: the diagonals of the caustics intersect at the 
same angle, since, as previously discussed, these angles are 
the same as 023,p- However the folds and cusps of these caus¬ 
tics look different. 

Having matched different lenses based on their 023,p 
value, we would like to see how the corresponding surfaces in 
the 3D angle space compare. We carry out the comparison in 
two steps: first we compare the edges and then the bodies of 
the surfaces. To carry out these comparisons we picked four 
different lens models, each belonging to the same 023,p = 92° 
contour. We identify each lens as L(e, 7 , P, 023,p). Table[^ 
summarizes these lenses. 

For the first comparison, the (180°, 034 , 023 ) edge of the 
lower of the two 3D surfaces for lenses A-D are plotted in 
Figure]^ Visually they seem to form the same line. But in 
order to quantify any difference between the edges, we com¬ 
pare fit equations of the lines. The fit equations for Lenses 
A-D are 

023 = 0.48617034 + 0.25271 

023 = 0.48617034 + 0.252017 , . 

023 = 0.48621034 + 0.247357 

023 = 0.48627034 + 0.240303 

respectively, and the angles are expressed in degrees. The av¬ 
erage of the median observational error of the three relative 
angles is ~ 0.8° (WW12) which is at least two orders of mag¬ 
nitude greater than the difference between the intercepts of 
the above equations. Similarly, the median observational er¬ 
ror of the slope 023/034 is ~ 0.006 which is again two orders 
of magnitude greater than the difference between the slopes 
of the four equations above. A corresponding comparison for 
the other edge of the 3D surface, the (0i2,180°, 023 ), yields 
similar results. Therefore, the difference between the edges 
of these four, and presumably other different SIEP-|-shear 
lenses in the 3D space is negligible for all practical and ob¬ 
servational purposes. 

For the second comparison we generate two sets of ran¬ 
dom source positions on one of the two surfaces of each lens. 
The quads from the first set are used to obtain an interpola¬ 
tion fit equation, which is then compared to the quads of the 
second set. This is our calibration, or cross-validation proce¬ 
dure. The comparison is done with the root mean square of 
errors, RMSE. The cross-validation for Lens A yields RMSE 
of 0.005°. The RMSE’s of Lenses B-D, when compared to 
the fit equation of Lens A are 0.011°, 0.029°, and 0.017°, 
respectively (see Table 1). A visual indication of how close 
these surfaces look is provided by Figure]^ which plots a 
certain projection of the 3D angles space (explained in de¬ 
tail in Section for Lenses A-D. The largest deviation is 
for Lens C, represented by (green) diamond symbols, which 
extend lower along the vertical axis by about 0.1 degree 
compared to the other three lenses. In general, lenses with 
largest deviations from the rest of their family members are 
those with larger ratios of e/ 7 . 

The RMSE values and Figure show that the surfaces 
are not mathematically Identical, but because the deviations 
between lenses sharing the same 023 ,p are smaller than the 


typical observational error (even for Lens C), observation- 
ally speaking these four SIEP-fshear lenses of different lens 
parameters (e, 7, P) are degenerate in the 3D angles space, 
l.e. in terms of the relative image angles. In other words, 
if a lens model described by (ei, 71, Pi) is found to fit the 
relative angles of an observed quad, then any other lens char¬ 
acterized by the same 023,p, but a different set (e2, 72, P 2 ) 
will also be able to reproduce the image angles, though for 
a different location of the source. 

So the answer to the question raised at the top of this 
section is yes, a single point in the 3D angles space, 023,p, 
adequately represents the rest of the surfaces. 

We note an interesting consequence of the degeneracy 
described above. Taking a (?=const two dimensional slice 
through the surfaces in Figure we get contours plotted 
in Figure These contours are symmetric about e = 7 
line, which implies that e and 7 are interchangeable. This is 
unexpected because this symmetry is not at all obvious from 
the equation of the lens potential, eq. with a = 1. Later 
we will see that this symmetry applies to other potentials if 
the isopotential contours are purely elliptical (as defined in 
Section 3.21. 

Note that in this paper we do not consider image dis¬ 
tances from the lens center; two lenses degenerate in image 
angles need not be degenerate in image distances. 


3.2 Non-isothermal Elliptical Potentials with 
External Shear 

As discussed when introducing the general power law po- 
tential, eq. [Tj a i s observationally constrained to be ~ 1 
(Koopmans et al. 20061. Therefore, we extend our discus¬ 


sion of Type H lenses only somewhat beyond this value; we 
use a = 0.9 and a = 1.2. We expand a up to 1.2 so that 
we can explore lenses with higher values of ellipticity. Going 
beyond these limits results in the mass density contours be¬ 
coming peanut shaped even for moderate values of ellipticity 
e. For 0 = 1.2 the two sheets in the 3D space starts crossing 
for small values of 023 when the ratio e /7 ^ 1 . 

We note that the way a is introduced in eq. creates 
lenses whose isopotential contours due to the main galaxy 
are not exactly elliptical; we call these hybrid potentials. 
The pure elliptical potentials are obtained by raising what 
we call the elliptical radius, rsjl — e cos 20 , to the power of 
a instead of raising just r. But the angular distribution of 
quads generated by pure elliptical potentials are not dis¬ 
tinctly different from that of SIEP-|-shear, therefore our dis¬ 
cussion in this section focuses on hybrid potentials. 

The 7 vs. £ contour is not symmetric about 7 = 6 line, 
with the contours rotated clockwise for a = 0.9 and counter 
clockwise for a = 1.2; see Figure]^ But for pure elliptical 
lenses defined above, the contours remain symmetric inde¬ 
pendent of a and this symmetry extends to NFW profiles 
discussed in the next section, as already alluded to at the 
end of Section 13.1.21 

We would like to know if the degeneracy property of 


SIEP-fshear (Section 3.1.2 1 applies to non-isothermal lenses, 
l.e. if these lenses can be adequately characterized by just 
023,p. For that purpose we generate quad distributions for 
Lenses F - I; see Table Lens F and G have a = 1.2, 
while H and I have a = 0.9 and all share the same peak of 
023,p = 93°. We use the second of the two tests Introduced 
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in Section |3.1.2[ We find fit equations for one of the two 
surfaces of Lenses F and H, and compared these to the actual 
surfaces of Lenses G and 1. 

The cross-validation RMSE for Lenses F and H are 
0.013° and 0.003°. The RMSE for Lens I as predicted by 
H is 0.264°. The RMSE for Lens G as predicted by F is 
0.551°. However, the high e to 7 ratio of Lens F (e /7 = 4.5) 
results in its two surfaces in 3D angles space crossing each 
other, which makes it hard to separate out a single surface 
and find a fit equation for it. Therefore the latter value can 
be compared only approximately to other RMSE values. 

This comparison of different lenses within the same fam¬ 
ily (F vs. G and H vs. I) shows that non-isothermal hybrid 
potential lenses sharing the same peak 623 ,p are approxi¬ 
mately degenerate in the 3D angles space (all RMSEs quoted 
above are smaller than observational uncertainty), but are 
not as similar to each other as those within the SIEP-|-shear 
family. 

We note, however, that this comparison might be af¬ 
fected by the fact that these calculations were done numeri¬ 
cally; SIEP-bshear potential is simple enough to be amenable 
to semi-analytical calculations, but non-isothermal hybrid 
potentials require fully numerical computations, and there 
are two places where numerical errors can creep in: (i) when 
selecting (e, 7 , /3) sets from the same 023,p contour, and 
(ii) when doing forward lensing of sources to produce im¬ 
ages and measure their polar angles. These numerical errors 
could contribute to the differences in the 3D surfaces. In the 
next section, where these errors are also an issue, we carry 
out a test to assess the error arising from (i). 

As with isothermal lenses, the models that differ the 
most from the rest of the family are those with larger e /7 
ratios. For a = 1.2 larger e /7 ratio also make the two sur¬ 
faces in the 3D angle space cross the FSQ and each other. 
Nevertheless, the differences quoted above are at most of the 
order of observational errors. Therefore observations based 
on the relative image angles of quads can not discriminate 
between any of the lenses within the family given by eq. 
as long as e /7 ^ 2 — 2.5. 


4 NON POWER LAW POTENTIALS AS 
EXAMPLES OF TYPE II 


In this section we broaden our exploration of Type II lenses 
to include models with non power law density profiles. We 
chose NFW radial density profile | | Navarro et al. I pffl), 
which has an analytical form for its lensing potential (Golse 


& Kneib 20021. Though the central density profiles of lensing 


galaxies are significantly affected by baryons and so typically 
have steeper than NFW slopes, we use this profile to explore 
how non power law profiles behave as Type II lenses. In or¬ 
der to make a Type II lens with NFW radial density profile 
we introduce external shear and ellipticity to the potential. 


<l>{r,0) = 61 In^ 


1 + VT^- 
1 - 


-f|r^cos(2[6i-/3])(3) 


where r = XyJ\ — ecos( 20 )/rs, x is the polar radius in the 
lens plane, is the scale factor of NFW density profile, 
and 6 is a normalization that is a function of and the 
characteristic density. The rest of the variables are as defined 


in previous sections. This potential is defined for r < 1, i.e 
for X < Xs- 

The diamond caustic of this potential has the symme¬ 
try of Type II lenses, and gives rise to two surfaces in the 
3D space of relative image angles of quads. As shown in 
upper right panel of Figure 7 vs. e contour map of con¬ 
stant 023,p is symmetric with respect to 7 = e line. Lenses J 
and K in Table are two different lenses in potential space 
whose two surfaces each have the same peak, 023 ,p, in the 
3D angles space. To check if these surfaces of two different 
lenses are degenerate in 3D angles space, the second test of 
Section [3.1.2| is implemented. A cross-validation RMSE for 
Lens J is 0.006° while the RMSE of Lens K as predicted by 
Lens J is 0.181°. 

Again, numerical errors could contribute to the differ¬ 
ence between the 3D angles space surfaces of these two 
lenses. To show that this is likely to be the case we made 
use of the symmetry with respect to the 7 = 6 line men¬ 
tioned in the previous paragraph. We generated 3D angle 
space surfaces for sets of lenses related by this symmetry: 
Lenses K and KK, and J and JJ in Table 1 are two exam¬ 
ples. All four lenses share the same 023,p = 93°, but the e 
and 7 parameters of J and K were picked from numerically 
generated output, whereas those of K and KK (J and JJ) 
were obtained by swapping e and 7 values, which yields ex¬ 
act e and 7 parameters. The RMSE from the comparison of 
Lenses J and JJ (K and KK) yields 0.007° (0.036°), values 
which are considerably smaller that 0.181° quoted in the 
previous paragraph. 

Whatever the source of the discrepancy between lenses 
such as J and K, the difference is still less than the observa¬ 
tional errors, so elliptical NFW-|-shear lenses with the same 
023 ,p are nearly degenerate in terms of image relative angles. 


5 IMPLICATIONS FOR THE OBSERVED 
QUADS AND SUBSTRUCTURE IN LENS 
GALAXIES 


In this section we rely on the basic property of Type II 
lenses, namely that the quads of each lens generate two dis¬ 
tinct surfaces in the 3D space of image relative angles, such 
as shown in Figure]^ (Section 3.1.11, and that the two sheets 
are qualitatively similar for all Type II lenses that are suit¬ 
able as models for observed quads. 


5.1 Why degeneracies are useful 

To draw conclusions about the broad class of lenses used as 
models one needs to explore a wide range of Type II lens 
models that are used in the literature to model quads. This 
task is made easier by the results of Sections |3]|^ where 
we explored detailed properties of three different represen¬ 
tative classes of parametric models as examples of Type II 
lenses, and showed that lenses within the same class, like 
SIEP-fshear, or elliptical NFW-fshear exhibit approximate 
degeneracies, in the sense that lenses with the same angle 
023 ,p have very similar distribution of quad relative angles 
in the 3D angle space. As mentioned earlier, we also verified 
that other commonly used profiles, most notably Sersic, also 
share these properties. 
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Furthermore, we now show that this degeneracy ex¬ 
tends to lenses with different radial density profiles shar¬ 
ing the same peak, 6 ^ 23 ,p, as long as the ratio e /7 is smaller 
than 2 — 2.5. To show this we introduce another lens. Lens 
E (SIEP-l-shear), which has the same 023,p as Lenses F 
- K. Interpolation fit for Lens E gives a cross-validation 
RMSE=0.003 degrees. Then, fitting Lenses G (hybrid with 
a = 1.2), H (hybrid with a = 0.9), and J (NFW) give RMSE 
of 0.146, 0.011, and 0.085 degrees, respectively (Table 1), 
which are all smaller or of the order of observational un¬ 
certainty. This is an important result: most Type II lenses 
which have the same angle between 2nd and 3rd arriving 
images of the central source are nearly degenerate in terms 
of their image angles, for all sources. 

The existence of this degeneracy, even if approximate, 
reduces the number of lens models one has to consider. It 
follows that we have to consider only a set of lenses with 
different 023 ,p, and any density profile (power law or curved 
in log-log space, like NFW, with the projected density slope 
not too different from 1 and any set of (e, 7 , /3), as long 
as e /7 ^ 2 — 2.5. The restrictions on a and e are because 
the degeneracies break down outside of the specified ranges. 
However, almost all of the lens models used to fit observed 
quad systems belong to the set of degenerate lenses. It is also 
important to stress that the general shape of the surfaces 
formed by quads in the 3D angles space is the same for all 
Type II lenses, even if they are not degenerate. 


(1) In order to reproduce quads at small O23 Type II 
lenses would need to have high values of ellipticity and shear. 
For example, the full SLAGS lens samples were fit with SIE 
mass profiles plus external shear model with median (max¬ 
imum) values of 7 and e of 0.05(0.27) and 0.21 (0.49), re¬ 
spectively ( Bolton et al.| 2008 1 , and the corresponding values 
from |Sluse et al. (20121 for 14 quads are 0.12 (0.33) and 0.20 
(0.49). The higher end of these parameter ranges appear to 
be larger than what should be expected in a realistic leasing 
case. Using monte carlo simulation and considering the ef¬ 
fects of environment Wong et al. (20111 have found the typi¬ 
cal value of tota l shear to be 0.08 with highest value of 0.17, 
and Dalai (20051 calculated a typical value for shear of 0.03|^ 
Similar results were obtained by Keeton et al. (19971 and 


Holder & Schechter (20031. These works imply that Type H 


lenses used as models for some observed quads have higher 
7 values than realistically expected. 

(2) In Figure the relative angles of all quads at 
|A023| < 12° and 023 < 45° can be modeled with Type 
H lenses shown in that plot as gray surfaces. However, these 
same Type H models also predict a large population of quads 
at 023 > 50° and | A023| > 5°. In fact, any Type H model pre¬ 
dicts that about 1/3 of its quads should have 023 >60°. This 
is in contradiction to the observations which show virtually 
no quads in that region. We conclude that even though Type 
H lenses are able to model most quads individually, they in¬ 
troduce a serious problem because they predict the existence 
of quads at large 023 and A 023 of at least a few degrees. 


5.2 Type II lenses cannot fit the population of 
observed quads 

Now we are ready to compare Type H lenses with the ob¬ 
served quad population. We will work in two dimensions 
instead of the three dimensions of the 3D angles space. We 
use the FSQ as a reference surface, and plot the vertical 
“distance” between a quad in the 3D relative angles space 
and the FSQ, called A023, versus 023. This plane is shown in 
Figure]^ together with the quads data (black filled symbols 
with error-bars) originally presented in WW12. FSQ in this 
plot is the A023 = 0 line. Using this same diagnostic tool. 
Type I lenses were shown to be inconsistent with the ob¬ 
served quad population (WW12). Specifically, the deviation 
of even extreme Type I lens models from the FSQ is not 
enough to account for the wide spread of observed quads in 
A 023 vs. 023 plane, especially at the lower values of 023 . 

Gray regions in Figure [^represent four examples (with 
two sheets per lens) of Type H lenses. Larger values of 7 or e 
move the gray surfaces further away from the FQS. Looking 
at the gray surfaces shown in that Figure, and remembering 
that many more similar surfaces can be drawn, both closer 
and further away from the FSQ, one may conclude that the 
whole wide class of Type H lenses—the most commonly used 
parametric lens models—should be able to reproduce the 
quad data. This impression comes about because one can 
find a Type II model, i.e. a set of (e, 7 , /3) that would 
reproduce any individual quad (with the possible exception 
of the extreme quads at |A023| ^ 15°), or in the parlance 
of Figure one can find a Type H model whose surfaces go 
through any given quad. 

However, this impression that Type Hs can reproduce 
the quad population is false, due to two reasons, with the 
second one being more important. 


5.3 Is our conclusion affected by a lens selection 
bias? 

Since the population of quads we use does not represent a 
homogeneous sample, one may wonder if selection biases are 
responsible for the distribution of observed quads in the 3D 
angle space and hence in the A023 vs. 623 plane. In particu¬ 
lar, can selection effects negate the argument in ( 2 ) above? 

To see that this is not the case let us recall how an¬ 
gle 023 of a quad is related to the source’s location within 
the caustic. Quads with large 623 originate from the cen¬ 
tral regions of caustics, while those with small O23, from the 
outer regions, adjacent to the folds and cusps. To account 
for the distribution of observed quads described in ( 2 ) one 
would need to have the central quads arising from Type H 
lenses suppressed by some selection bias, while not suppress¬ 
ing the central quads arising from Type I lenses. It would be 
hard for a selection bias to accomplish that. For example, 
magnification bias would bias would suppress central quads 
because they have smaller magnifications. However, it could 
not distinguish between Type I and Type H lenses, which 
is what would be needed to explain the distribution seen in 
Figure 

5.4 Modeling individual quads vs. modeling quad 
population 

We stress that our analysis presented above differs signifi¬ 
cantly from lens modeling of individual quads. Unlike lens 
mass modeling, parametric or free-form, which is done with 

^ Note that these ellipticity values refer to mass, so cannot be 
compared directly to our e. 
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one lens system at a time, we study the population proper¬ 
ties of quads, and compare them to the generic properties 
of Type II lenses. This is why we refer to our method as 
model-free analysis of quads. As pointed out above, and can 
be seen from lens modeling literature, most individual quads 
can be successfully modeled with Type II lenses. However, 
these same lens models also predict the existence of quads 
which are not observed, and whose absence would be hard 
to explain by a selection effect. This predicted, but unob¬ 
served population of quads is most clearly seen in the space 
of relative images angles, as is done in Figure 

5.5 Type III lenses 

This leaves us with Type III lenses, for example substruc¬ 
tured lenses. A preliminary analysis gives promising results. 
As depicted in Figure]^ a better match to the distribution 
of quads in the A623 vs. 623 plane is achieved with the intro¬ 
duction of substructure in the main lens. By increasing the 
dumpiness of the lens, the gaps in the peak region, near large 
023, of the two surfaces of quads decreased significantly and 
looks more like the observational data. However, too much 
dumpiness could lead to higher order catastrophes in the 
caustic which would result in observationally unsupported 
systems with higher image multiplicities. In the examples 
used in Figure this is not a concern since even the clump¬ 
iest of the two substructured lenses (bottom row) produce 
more than five images for less than 1 % of the total sources. 

Figure [T0| shows that one does not need to use a highly 
substructured lens to reproduce the population distribution 
of observed quads. Here we start with a substructure-less 
Type I lens (top row) that gives rise to quads which are 
closely packed on the horizontal axis, i.e. the FSQ, failing 
to match the spread of observed quads. The introduction 
of substructure (bottom left) does not perturb the shape of 
the outer isodensity contours, but the resulting distribution 
of quads is dispersed in the A023 vs. 623 plane in the same 
fashion as the observed quads (bottom right). 

This is just a limited foray into Type HI lenses. The 
conclusion of this Section is that Type H lenses cannot re¬ 
produce the distribution of image angles of observed quads, 
and that possibly substuctured lenses, or additional nearby 
perturber galaxies are needed. We leave quantitative work 
on these lenses, and detained comparison with the observed 
quad distribution for a later paper. 


6 CONCLUSIONS 

In this paper we presented a technique to explore the prop¬ 
erties of quadruple image gravitational lens systems. In gen¬ 
eral, the relative distribution of quad’s four images about 
the center of a lens can be fully described in 6 D space of 
three relative angles and three relative radial distances in 
the lens plane. Our technique uses only the 3 D subspace 
of relative angles which is orthogonal to the remaining 3 D 
subspace of relative distances. 

We classify lenses into three main categories based on 
whether their diamond caustic obey twofold and double mir¬ 
ror symmetries. The caustics of Type I lenses fulfills both 
symmetries, those of Type H obey only the twofold sym¬ 
metries while the caustics of Type His obey neither. In this 
paper we focused on Type H lenses. Our main working space 


is the 3 D space spanned by the three relative image angles 
of a quad, 612, O23, and 634,. One of our main conclusions is 
that this space is a useful tool in studying quad lenses. We 
showed how the distribution of the three angles of a given 
lens relates to its leasing potential and the caustic. The two 
fold symmetry of the Type H lens caustic is reflected in the 
fact that in the 3 D angle space the distribution of the rel¬ 
ative angles is always confined to two surfaces (Figure]^. 
This is the defining feature of Type H lenses. 

Of the three potentials we studied the closest quanti¬ 
tative connections between the lens potential, caustic and 
3 D angle space are exhibited by the SIEP-|-shear potential 
(eq. [^with a = 1 ). The angle between the second and third 
arriving images of the source located at the center of the 
lens, 023,p is equal to the angle between the caustic diago¬ 
nals. (Though we do not study the images distances from 
lens center in this paper, we note that the distance ratio 
of these two images is equal to the ratio of the two caus¬ 
tic diagonals.) SIEP-|-shear lenses that have the same angle 
023,p have nearly the same distributions of quads in the 3 D 
angle space, implying that 023,p is sufficient to specify the 
shape of the surfaces of all these lenses, and that these sur¬ 
faces are degenerate. Other lenses also show close similarity 
between lenses of the same 023 ,p, but not as close as for 
SIEP-fshear. The similarity extended to lenses of different 
radial density profiles: quads of Type H lenses with elliptic- 
ities e/7 < 2 — 2.5 that have the same 023,p share approxi¬ 
mately the same two surfaces in the 3 D angles space. 

In addition to the lens potentials discussed in detail in 
this paper we explored more general forms such as the Ser- 
sic profiles. The commonalities in the distribution of quads 
in the 3 D angle space (the general shape and location of 
the two surfaces) persist in all the models considered. This 
should not come as a surprise since it is a direct reflection 
of the two-fold symmetry of the caustics of Type Hs which 
identify the class. Furthermore, the near degeneracies de¬ 
scribed above persist across all the models studied and we 
conjuncture that they extend to all Type Hs since there is no 
evidence to the contrary. Because these observations apply 
to a wide range of models and allow us to draw conclusions 
about Type H lenses as a class, we call this approach model- 
free analysis. We conclude that, even though Type H lenses 
can successfully model individual quad lenses, they cannot 
reproduce the relative image angle distribution of the pop¬ 
ulation of the observed quads. 

We present two examples of Type HI substructured 
lenses that have distributions of relative image angles that 
are generally consistent with the observed distribution. 
These examples are not meant to be a close match to obser¬ 
vations. Future work will need to look at different types of 
substructure, as well as luminous and dark secondary per¬ 
turber galaxies more carefully to ascertain what type and 
amount of substructure is necessary to fit the observed quad 
population. Finally, we note that the substructure mass 
clumps one would need to explain image angles are likely 
to be much larger in mass and extent compared to those 
thought to be responsible for the flux ratio anomalies of 
quads. 
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Figure 1. Lens Parametrization in three types of spaces: Left: Potential space. A Type II elliptical lens potential contours with external 
shear oriented along the line (inclined, and going through the center). The four images are numbered based on their arrival time Middle: 
Caustic space, parametrized by the ratio of its diagonals (red lines) and the angle they make at their crossing. S is the source that gives 
rise to the images shown in the left panel. Right: 3D space of quad relative image angles: three dimensional space formed by the three 
relative angles 9 i 2 , O 34 and O 23 . Each point in this space represents a quad. The surfaces are generated by images of many randomly 
distributed sources within the diamond caustic and the cut. 



Figure 2. 3D angles space: Type II lenses produce quads that form two sheets on either side of the FSQ, shown as the smooth surface in 
both panels. (Note that it looks dark red from the left and bright red from the right). Quads belonging to the Type II lens are represented 
by light blue dots. Left: The top portion, 623 > 75°, of the 3D space and the surfaces. The peak of the FSQ is at 023,p = 90°. The peaks 
of the two surfaces of Type II lens are at 023,p ~ 94°, and 023,p ~ 86°. Right: The lower portion of the same 3D space, for 623 < 20°. 
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Figure 3. Mapping of the caustic space to the 3D space of relative image angles. Top: Each pair of opposite quadrants of the diamond 
caustic (for example, the two black quadrants) map into one of the two surfaces in the 3D space (black). This is a reflection of the twofold 
symmetry of Type II lenses. Bottom: Sources along constant position angle in the caustic plane (various colored lines) are mapped to 
non-crossing curves in the 3D angles space (similarly colored lines). 



Figure 4. Contour surfaces of constant ^ 23 , p for SIEP-|-shear lens in 7 vs. e vs. j3 space. ^ is in radians. The outer (inner) surface is for 
023,p = 91° (^ 23 ,p = 92°). The jagged line corresponding to ^ = 0 and 7 = e is the result of no quads being formed for those parameters. 
Two orientations of the same space are shown. 
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Figure 5. The 612 = 180° edge of the lower of the two surfaces, in the 3D space of relative angles, for four SIEP+shear lenses sharing 
the same peak 023,p- Left: Bottom end of the edge. Right: Top end of the edge. The parameters of the four lenses A-D are given in Table 
1 . 
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Figure 6 . Pojection of Lenses A-D with parameters as given in Table [T], Lens C, diamond (green), slightly deviate from the other three 
lenses by saging at the folds. 


© 2015 RAS, MNRAS OOO.plfTsl 




14 Woldesenbet and Williams 


SIEP+ Shear NFW 



Figure 7. Contours of constant 023,p in 7 vs. e plane for lenses of four different radial density profiles. The contours are labeled by 
023,p- For SIEP+shear (Top left) and NFW {Top right) lenses the contour lines are symmetric about 7 = e line. For hybrid power law 
potentials the contours are rotated, as compared to that of SIEP+shear, clockwise for a < 1 {Bottom left) and counterclockwise for 
Q: > 1 {Bottom right). 
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Figure 8 . Distribution of observed quads (black filled circles) and quads from four Type II lenses (gray distributions) in 623 vs. A 623 
plane. The quad data and the errorbars were taken from WW12. A023 = 0 is FSQ. Each lens has two gray surfaces which are about 
equally close to, but on opposite sides of the FSQ. From the outermost to the innermost surfaces, the lens are: NFW (0 = O.lrad, 
e = 0.27,7 = 0.25), Power Law (a = 1-2,/5 = 0.2rad, e = 0.25,7 = 0.25), Power Law (o; = 1 ,^ = 0.2rad, e = 0.15,7 = 0.17), and Power 
Law (a = 0.9, 0 = 0.15rad, e = 0.14, 7 = 0.13), respectively. 
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Figure 9. Left panels: All three have isothermal radial density profile with ellipticity e = 0.25 and shear 7 = 0.2 oriented at an angle 
= 7 r /6 to each other. The k = 1 contour is shown by the dashed, red, contour. Right panels: The corresponding distribution of quads 
from these lenses in the 623 vs. A023 plane are shown as the the light shade, yellow, distributions. The observed quads (same as in 
Figure]^ are red (black) dots. Top: A Type II lens with no substructure. Middle: A Type III lens with randomly distributed clumps 
that account for 4.38% of the total mass within the window. Bottom: A Type III lens with randomly distributed clumps that account 
for 17.96% of the total mass within the window. 
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Woldesenhet and Williams 


Lens 

name 

Radial Profile 

L (/J, e,7i®23,p) 

Cross-validation 
RMSE (deg.) 

Lens comparison 
RMSE (deg.) 

A 

SIEP 

L (1.2, 0.09, 0.10366, 92) 

0.005 


B 


L (1.2, 0.06603, 0.19, 92) 


A: 0.011 

C 


L (0.5, 0.14, 0.036752, 92) 


A; 0.029 

D 


L (0.5, 0.034944, 0.06, 92) 


A: 0.017 

E 


L (f, 0.069197, 0.08, 93) 

0.003 


F 

Power Law 0 = 1.2 

L (|, 0.2, 0.0446, 93) 

0.013 


G 


L (J, 0.06536, 0.18, 93) 


E: 0.146, F: 0.551 

H 

Power Law 0 = 0.9 

L (|, 0.04706, 0.11, 93) 

0.003 

E: 0.011 

I 


L (|, 0.095, 0.058, 93) 


H: 0.264 

J 

NEW 

L (|, 0.054, 0.075, 93) 

0.006 

E: 0.085 

JJ 


L (|, 0.075, 0.054, 93) 


J: 0.007 

K 


L (|, 0.13, 0.05265, 93) 

0.002 

J: 0.181 

KK 


L (|, 0.05265, 0.13, 93) 


K: 0.036 


Table 1. Examples of Type II lenses. L (/3, e, 7 , 623 p) refers to a lens with ellipticity e, external shear 7 , which is oriented at an angle of 
/3 (in radians) relative to the principal axis of the main lens, and relative image angle 023,p between images 2 and 3 of a central source. 
The subscript p in 023,p indicates that it is the 623 value of the quad located at the peak of the top surface in the 3D space of relative 
angles. The cross-validation column lists RMSE values between an interpolation fit to the lens and its own quads. The lens comparison 
column lists RMSE value between the fit to a lens (listed in column 5) and quads from a different lens (lens named in that row). Lenses 
JJ and J (and also K and KK) have their e and 7 values swapped. 

This paper has been typeset from a T^K/ DTgX file prepared by the author. 
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